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Abstract 

It is well known that the correlation functions of a scalar field in a quasi-de Sitter space 
exhibit at the loop level cumulative infra-red effects proportional to the total number of 
e-foldings of inflation. Using the in-in formalism, we explore the behavior of these infra- 
red effects in the large limit of an 0{N) invariant scalar field theory with quartic self- 
interactions. By resumming all higher-order loop diagrams non-perturbatively, we show 
that the connected four-point correlation function, which is a signal of non-Gaussianity, is 
non-perturbatively enhanced with respect to its tree-level value. 
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1 Introduction 



In the inflationary picture the primordial cosmological perturbations are created from quan- 
tum fluctuations which are "redshifted" out of the horizon during an early period of accel- 
erated expansion [1]. These super-horizon perturbations remain "frozen" (unless isocurva- 
ture fluctuations are present) until the horizon grows during a later matter- or radiation- 
dominated era. After falling back inside the horizon they are communicated to the primor- 
dial plasma and hence arc directly observable as temperature anisotropies in the Cosmic 
Microwave Background (CMB) radiation. 

Given the spectacular accuracy by which these anisotropies are and will be mapped by 
present and forthcoming satellite experiments, the issue of non-Gaussianity (NG) in the cos- 
mological perturbations has becoming more and more relevant [2]. Such non-Gaussianities 
are sourced by self-interactions, including the ones induced by gravity, in the early universe, 
and become manifest at the level of second- or higher-order perturbation theory. Further- 
more, and more important for the scope of this paper, self-interactions of any scalar field 
during the inflationary stage give rise to corrections in all correlators, including the ob- 
servationally interesting cases of the two- and the three-point correlation functions. Such 
corrections are associated with so-called "loops," in which virtual particles with arbitrary 
momentum are emitted and re-absorbed by the fields which participate in the correlation 
function. 

While loop corrections may lead to significant backreaction effects [3-7] which scale like 
(powers of) the number of e-folds between horizon exit of the mode k under consideration 
and the end of inflation, in this paper we will concern ourselves with another aspect of these 
loop corrections. Indeed, it has been known for a long time that there exists a breakdown 
in the perturbative expansion due to infra-red (IR) "divergences" [8,9]. These classical one- 
loop corrections scale like ln(/cL) where is the IR comoving momentum cut-off. These 
IR divergences have attracted lot of attention in the literature recently, where in particular 
the loop corrections to cosmological correlation functions has been studied [10-16]. At 
the generic n-th order of perturbation theory, the power spectrum of the perturbations is 
expected to get corrections growing like ln"(A;L), signalling the breakdown of perturbation 
theory unless the corresponding coupling constant is tiny enough. 

One option is to set L to be some comoving length scale which left the horizon many 
e-folds before the observable universe (super- large box). In particular, the smallest possible 
value of is aiH, where Oj is the value of the scale factor at the beginning of inflation. 
Since the wavelength k^^ goes out of the horizon when the scale factor equals = k/H, 
we see that ln(/cL) = ln(afe/aj), which is proportional to the total number of e-folds from 
the beginning of inflation to the time when the mode k exits the horizon. This is the large 
cumulative IR correction as it gets larger the longer inflation lasts. Indeed it was estimated 
that in models of chaotic inflation, where the total number of e-folding is large, the one-loop 
correction to the power spectrum can be significant [10,11,14-16]. This may have important 
indirect implications for our understanding of the relation between observational quantities 
such as the tensor-to-scalor relation and the underlying fundamental model [10,11,14-17,20]. 
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If we insist in adopting a superlarge box to cut-ofF the IR divergences, it should be kept 
in mind however that the location of the much smaller box which includes our observable 
universe may be untypical and one should quantify how likely it is that the correlators 
averaged over the superlarge box coincide with the correlators in our observable universe 
[19,20]. To deal with this problem, one can try to use the approach of stochatic inflation 
[18] where the uncertainty in the prediction coming from the large IR contributions in the 
superlarge box is traded for the uncertainty inherent in having a probability distribution for 
the background quantities [20]. However, before attacking the problem of the IR divergences 
from this point of view, we should first ask ourselves whether it is possible to deal with 
them using the traditional quantum field theory techniques. We may try to resum them, 
for example by using the standard method of the Wilsonian Renormalization Group (RG). 
Unfortunately, this technique is not very efficient in practice. Indeed, the full set of exact 
RG equations need to be solved. Approximations inevitably lead to disregarding a set of 
diagrams which are not at all subleading in the IR. We refer the reader to the Appendix of 
Ref. [20] for more details. 

We are though aware of one quantum field theory example in which the resummation of 
the IR divergences may be performed, that is the case of N scalar fields subject to an 0{N) 
symmetry and with a quartic interaction in which the free parameter is taken to be large. 
This toy-model has been proved to be very useful in dealing with similar problems of IR 
resummation of perturbation series in the case of finite temperarure field theory [21,22] and 
non-equilibrium field theory [23,24]. We will then review how in the large A^ limit the 0{N) 
theory embedded in de Sitter can be solved and that the IR divergences can be resummed. 
Through the in-in formalism [25], one is able to give the expression for the exact two-point 
correlation function of the perturbations of the A^ scalar fields and to compute the exact 
connected four-point correlation function, dubbed trispectrum, which signals the presence of 
NG. The resummed trispectrum will show a parametric dependence on the quartic coupling 
which differs from the tree-level one. More important, the correlation functions will not 
manifest IR divergences. Despite the fact that we will be dealing with a toy model, the 
results obtained in this paper may be useful to understand the behaviour of more realistic 
systems in the IR. 

The resummation of IR divergences was previously discussed by Hu and O' Connor in the 
2PI (Two-Particle-Irreducible) formalism [26], which we review in section 2 and 3, and by 
Starobinsky and Yokoyama using stochastic field theory methods [27] , the results of which 
we also review in section 3. 

Thus, the paper is organized as follows. In Section 2 we review the details of the model 
and the in-in formalism. In Section 3 wc review how to solve the toy model in the large 
A^ limit. In Section 4 we present the nonperturbativc computation of the trispectrum. In 
Section 5 we compute the correction to the comoving curvature perturbation coupled to A" 
scalar fields, which is the toy model studied by Weinberg in Ref. [13]. Finally, our conclusions 
are given in Section 6. 
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2 The 0{N) invariant model in the in-in formahsm 

We consider a 0{N) invariant generalization of a scalar field theory with quadratic self- 
interaction of the type X(f)^. The generalized 0{N) invariant field theory consist of scalar 
fields (f)i, i = 1,...,N, with an action which is invariant under the A^-dimensional real 
orthogonal group 

^=1 (d.cPid^cP, - m^cP.cP, - -^{cpicpif^ , (1) 

where summation over repeated indices is assumed. We will later embed this theory in an 
expanding universe. In the in-in formalism, the generating functional, Z, can be written 



Z[J+, J_,p(ti„)] 



(2) 



Yl j P</)i+P</)r exp i j dt' j d?xyf^{C[4)i+] - £[</.r] + Ji+4^i^ - Ji-cpi') 

j J ctp \. Jtin J 



The closed time path (CTP) is defined by the condition = <p^{t) and the boundary 

conditions in the asymptotic infinite past. Now, by the usual Gaussian integration, one can 
write the free field generating functional in terms of the propagators 



such that for the free field generating functional, one then has 



Zq oc exp 



^ 9 jfCfij Jj 



where Jf = {J^ , ) and the Green functions are given by 

Goi;{x,y)=i{Ux)Uy))^'\ 
G4r^x,y)^i{ct>{y)^ct>{x),)^'\ 



G'o++(x, y) = i{Tct>{x)ict>{y)jr^ = e{xo - yo)Goi^(x, y) + 9{y^ - x^)Golr{x, y), 
Go~{x, y) = i(t(/)i(x)0,(y))(°) = e{xo - yo)Gotj''{x, y) + ^(yo - xq)Go-Mx, y), 



(3) 



(4) 



(5) 
(6) 
(7) 
(8) 



where the sub or superscript (0) denotes the free field correlation functions. 



2.1 One-loop logarithmic divergence 

As mentioned in the introduction, it has previously been shown that the logarithmic IR 
divergences to the power spectrum of inflationary perturbations can lead to large effects, 
which ultimately requires a non-perturbative treatment [10,11,14-16]. Let us first briefly 
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review how these IR divergences manifest themselves in the present case. The perturbative 
expression for the two-point function is 

J ctp 

We can then calculate the one-loop contribution by expanding the exponential to first order 
and perform the required Wick contractions 

{MVo,k)UVo,k')) = i2nf6^'\k + k') j -J-^2Q [G,^(k- ,r^o,v)Go^ik- ,r^o,v)] 

x\Gq+{x,x) . (10) 

On the right hand side we have suppressed the latin indices on the Green functions. The 
expression above has two logarithmic divergences. The most severe comes from 

-j-pm. (11) 

where the power spectrum of the scalar fields, V^{k), is scale independent on super-horizon 
scales in de Sitter space. The integral therefore diverges logarithmically with the lower limit 
of the k integral. Naturally one may take the lower limit of the integral to be given by 
the largest scale which can fit inside the in inflationary patch at the beginning of inflation 
k = ttiHi. This implies that the logarithmic divergency will scale like the physical size of the 
inflationary patch or equivalently like the total number of e-folding of inflation. The other 
logarithmic divergence comes from the time integral in Eq. flTUl) . It depends on krjQ and scales 
like the number of e-foldings left of inflation after the mode k left the horizon. This time- 
dependence simply reflect that the perturbation variable is not conserved (the associated 
comoving curvature perturbation in the case of only adiabatic modes will be conserved, see 
Refs. [14,15]). Thus, the exact non-perturbative two-point function is expected to depend 
on • In order to deal with the former more severe IR divergency non-perturbatively, we 
will apply the Two-Particle-Irreducible formalism. 



2.2 Two-Particle-Irreducible formalism 

From the generating functional defined above, we can define the generating functional of 
connected Green functions W = —i log Z, and by a Legendre transformation, we can then 
compute the effective action F, which is the generating funtional of One- Particle-Irreducible 
(IPI) Green functions. 

However, in the non-perturbative regime where quantum fluctuations are large, the IPI 
effective action is insufficient because it only gives the dynamics of the mean field 0c = {4>) 
and does not describe the evolution of the variance (0^). In the 2PI approach the functional 
derivative of the effective action by the Green functions yields an equation for the variance. 
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called the gap equation, just like the dynamics of the classical field is given by the functional 
derivative of the effective action by the classical mean field. 

In the 2PI approach, one generalizes the generating function to include a non-local source 
term [28, 29] 



n/ 

i,a "^"'P 



L>0iexp 



where summation over repeated indices is implicit, the matrix Cab is defined as c 
1, c+~ = = 0, and 

Sm= [ riS;v/=^(/:[0>]-£[0r]). 



++ 



(12) 



-c = 



(13) 



The effective action is a double Legendre transform of W{J,K). One defines the classical 
field and the two-point function as 



"^SJ^ix 



cf>,t{xWj{x') + Gt}{x,x')^2 



6W 



-g^dK'^i^x.x') 



(14) 



(15) 



Now one can eliminate and K^j in terms of and G^j*, and define the 2P1 effective 
action as a double Legendre transform of W, 

r[<p,,G] = W[J,K] - J d^X^CabJt(t>c]5'' 

-If d'x^ f d'x'^CabCcdKt^{x,x')\Gl'i{x,x') + (l>Jl{x)^^^^^^ 



In a Feynman series expansion of F, one then obtains 



F[0„G] = 5[0e]-^lndet[Gg'] 



+ 2 / d'^^ I d^x'^G-,ZW.x)G'${x,x') + F^I^cG], (17) 



where 



5^S 



-g5ct>1{x)mx')^- 



(18) 
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After shifting the field around the background field (p ^ tjj in the classical action S[(j)] order 
to define a fiuctuation field ip = (f) — (f)c, the new action becomes S[(f)c + ip], and T2 is then 
given by all closed 2PI graphs (all 2-point graphs which does not become open by opening 
two propagator lines), with all the propagator lines given by G and vertices given by the 
shifted action. 

The mean field equation 



and the gap equation 



can then be calculated to yield 



Cc — f c =Vc 



SGt^{x,y) 







0, 



(19) 



(20) 



"c — Vc =Vc 



+ + 4 mx)+Gt;{x,x)]^ <t>.,{x) = ^<P.,ix)Gj^M - ^g^, (21) 



and 



Gt^M = ^4>c^{x)4>c,{x)Gt^{xy) (22) 
+ i j d^z^++{x,z-Ae,G)G+;{z,x') + t6,/-^^^0^ 



with 



^ti [x, y, <f>c, G) = 2z— p^- -. 

6Gj^{x,y) 



(23) 



3 Solutions in the large N expansion 

The leading order large approximation corresponds to a self-consistent mean field approx- 
imation, where the infinite hierarchy of Schwinger-Dyson equations for the n-point functions 
are truncated to a closed system of two dynamical equations of just the one-point function 
{(f)i{x)) and the two-point function {(j)i{x)(l)j{ii)) . Since no irreducible higher order correlation 
functions appear in the large N limit, it is equivalent to a Gaussian approximation for the 
density matrix. 

In order to calculate the leading limit of the mean field and gap equations, we note 
that a trace over the 0{N) indices gives a factor of A^, while the vertices each contribute a 
factor 1/A". One can for simpficity rescale the field and composite field operators [28,29] 

<l>,^{x) ^ VNrdx), (24) 
G^^{x,x') ^ G''\x,x%j, (25) 
ij^{x) ip'^ix) (26) 
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in order to obtain the mean field equation to leading order in the large N expansion 

-(n, + m' + ^ [(jylix) + G++(a;,a;)]^ M^) = 0. (27) 
Using that to leading order in large N 

^^""[^^ ^ "8^ / "^'"^^ [GS+(x,x)G++(a:,x) - Gr.-(a;, x)G,T(x, x)] 5^'^ (28) 
one can also obtain the gap equation for the rescaled fluctuation Green function 

-(D, + m' + ^ [02 (x) + G++(a;, x)] ^ G++(a;, x') = z '^^^'^ (29) 

In the equation above, we have omitted the appropriate counter terms needed to cancel the 
UV divergences, since in the following we will only concern ourselves with the IR properties 
of the equation. 



3.1 Large wavelength approximation 

The flrst step in solving the gap equation for G^~^{x,x') is to solve for G~^~^{x,x). Once we 
have an analytical expression for G~^~^{x, x), we can insert it in the equation above and solve 
for G++(a;, x'). We now assume that the N scalar fields are evolving during a de Sitter stage 
characterized by a constant Hubble rate H. In the large wavelength approximation, we can 
approximate the d'Alembertian by SHdt- Assuming further that initially (pc = 0, and that 
it will stay small due to the — > —0 symmetry of the Lagrangian, in the massless case the 
gap equation will at late times simplify to 

- 3HdtG++{x, x) = X [G++{x, x)] ' + ^2-^ , (30) 

where the extra factor of two is a symmetry factor, which accounts for the difference in taking 
the derivative of G{x, x') before or after taking the coincidence limit. The appearance of 
the ambiguous 5(0) implies that, in order to solve this equation, we must supplement the 
equation with a boundary condition. In the Hartree approximation in Ref. [27], it was 
obtained by requiring that one recovers the correct linear growth of the variance in the limit 
. 

One may however also observe, that the contribution from the IR modes is captured by 
introducing a smearing with a characteristic size corresponding to a fixed physical distance 
of the order of the horizon scale, and this gives rise to an additional term in the equation 
above governing the time evolution of the smeared two point function, as was also discussed 
in Ref. [27]. Below, we will instead solve for G++(x, x) using the Fokker-Planck equation 
and then stick the result back into the original gap equation and then subsequently solve for 
G++{x,x'). 
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3.2 Fokker-Planck equation 

The stochastic approach embodies the idea that the IR part of the scalar field may be 
considered as a classical space-dependent stochastic field satisfying a local Langevin-like 
equation. The stochastic noise terms arise from the quantum fluctuations which become 
classical at horizon crossing and then contribute to the background [27]. The expectation 
value of any function F[(j)\ of the coarse-grained field is determined by 

(F[0]) = j d<PF{<P)P{<P,t). (31) 

Using that in the massless case, the classical Lagrangian corresponding to the rescaled field 
in eq. is 

'C=^f9^09>- VI , (32) 



2 " 4 

one can derive a Fokker-Planck equation describing how the probability of scalar field values 
at a given spatial point evolves with time 

Now, let us consider the expectation value of the 2?7.th power of the field. Following Ref. [31], 
one can differentiate the expectation value of c/)^"' and use the Fokker-Planck equation to 
obtain 

I ¥'") = """4;/"'' "-^> - ^ (*"*^> . (34) 

if one assumes that the field vanishes at the boundaries of the field integration. Like in 
Ref. [31] one can define 

a = In a, A = An'^X/S, (35) 



47r2 

in order to obtain the differential recursion relation 

It can be solved iteratively to yield the perturbative solution of Ref. [31] 

n , . . A 



(0-> = (2n-l)!!(glna 



1 - -(n + 1)— ^ In^ a + ... 

2^ ^127r2 



(37) 



However, in the case we are interested in we can solve the equation exact. Using that in the 
large limit (0^) = (0^)^ (only the trace squared part contributes in large A^), we obtain 
for n = 1 

' - • (38) 
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which has the solution 




Tanh(^lna) 



(39) 



This agrees with the asymptotic solution for Ina oo given in Ref. [27], and exemplifies 
how the logarithmic IR divergences can be resummed to yield a finite result. The solution 
for (0^) = G^^{x,x) can now be inserted back into the gap equation, which can now be 
solved for G^^{x,x'). 



3.3 Analytical solution to the gap equation 



Let us assume that G^^{x,x) is given by Eq. ( l39i) and insert that into the gap equation 
Eq. ( l29l) and solve for G{x, x')~^~^ in the large wavelength limit, with x ^ x' . Then we obtain 
an equation of the form 



dt + ^ 



X tanh(at) 



G{t) = 0, 



(40) 



where a = y/XH/Air"^ and b = 3^/\/H. The solution is 

G(t) = v4(cosh(at))-^/'"V 



(41) 

The solution scales like exp[— (\/3A/127r)ift] at large times. If we are only interested in the 
asymptotic late time behavior, there is a more illuminating way to proceed. At late times 
the variance (0^) — > H'^/VJ^, and in this case the gap equation in Fourier space yields 



(42) 



+ 3Hdt + ^ + G'++(t, to; k) = 0, 



47r 



where we neglected only the gradient. Now the interpretation is clear. The non-perturbative 
effect is a non-perturbative regulating mass given by m^p = ^^^H^. We can solve the 
equation by the usual method of going to conformal time and rewrite the equation on the 
form 
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G++{r],W,k) = 0, 



where 



m 



np 



(43) 



(44) 



We defined the rescaled Green function as G = G/a. As a function of rj the two independent 
solution scales like —t^h'^P [—kr]) and y/^rjHl^\—kr]) . Thus the asymptotic late time be- 
havior scales as (— fc?])^/^"'' at late times. Thus, instead of going to a constant, G~^^{r], r]o; k) 
is decaying at late times as (— fc?^)™np/^i' the limit — > 0. This is the same late time 
behavior as given by the solution in Eq. fHTj) . The exponential decay in physical time of 
the retarded Green function, implies that the system has a finite correlation time, which is 
proportional to 1/VXH. 
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4 Non-perturbative enhancement of non-Gaussianities 



Let us consider the four point function of the 0{N) field in Fourier space. The interaction 
vertex for the rescaled fluctuation field in eq. ( 12^ is {X/8N)ip'^. At tree-level, the four-point 
function is given by 

(45) 

where Go is the Green function of the free field. The Green function of a free massless field 
in de Sitter is given by 

Go +(77, r/'; k) = t^ikv - i){kv' + t)e''^^'-^^ . (46) 

Expanding the real and complex part of the free field Green function to leading order in the 
large wavelength limit, one then obtains a expression for the four-point function, which is 
logarithmically divergent at late times. At late times the tree-level four-point function is 
(we simplified the expression isolating the leading terms) [32] 

{A, ...A.)^ ^^6i">{h + ■ • • + ^4)Jf^^^^^i^(7 + log(-fc,)). (47) 

The divergent log term is proportional to the number of e-folding left before the end of 
inflation after the mode has crossed the horizon, and for observationally interesting modes, 
it will be of order oneQ. Like for the two-point correlation function, one expects cumulative 
effects at one-loop coming from the IR. It is interesting to know how this conclusion is altered 
when we take into account all of the resummed loop effects by using the non-perturbative 
slotution for the Green function. 

The effect of the non-perturbative Green function is to regulate the log-divergent integral 
by the non-perturbative mass Tn"^^. With the non-perturbative mass, the Green function at 
late times take the form 

G{r], r]'; k) = t^{7]7]')^/^H^Hl^\-kr])Hl^\-kr]'). (48) 

If we expand the Hankel functions around zero, we can straightforwardly calculate the four- 
point function at late times. However, a less tedious way to proceed is to make an appropriate 
field redefinition. 

*At any rate, one has to take into account that, when computing the trispectrum of the gauge-invariant 
isocurvature perturbations, the time dependence may differ from the one in Eq. (|47p . 
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4.1 Change of basis and classical approximation 

Since we are interested in investigating the infrared properties of the inflaton propagator 
on super Hubble scales, where the dynamics can be described classically, it is convenient to 
make a rotation in field space, in order to make the comparison with classical field theory 
simpler. The change of basis is defined 



^(1) 

^(2) 



R 



R 



1/2 1/2 
1 -1 



(49) 



In the new basis the Lagrangian density >C ['?/''''] ~ ^Clip ], with the background field put to 
zero for simplicity, becomes 



^ ( _ ^2^(1)^(2) _ A (4(^(1))3^(2) ^ ^(1)(^(2))3) 

8N 



/:[^a),^(2)] = 

and the propagator matrix becomes 

GK(x,y) = RG(x,y)RT 

with 



iF{x,y) G^{x,y) 
G^{x,y) 



(50) 
(51) 



F{x,y) = --{G-+{x,y) + G+-{x,y)) , (52) 

G^(x,y) = G++{x,y)-G^-{x,y) = 9ixo-yo)iG-+ix,y)-G^-ix,y)) , (53) 
G^(x,t/) = G++ix,y)-G-+ix,y) = 9{yo-Xo){G^-{x,y)-G-+ix,y)) , (54) 

(55) 

Like in Ref. [33], we represent the -ip^^^ field with a full line and the tp^"^^ field with a dashed 
line. The Feynman rules for the two point functions and the vertices are then 



n T2_ 

n r2_ 

T2, 



T3 



1^2 



F{k,Ti,T2), 

-tG''{k,n,T2) = -iG^{k,r2,n), 
3A 

-i— a^(ri)5(ri - T2)6{ri - r3)5(ri - n), 
3A 



(56) 
(57) 
(58) 
(59) 



^For alternative diagrammatic approaches, see Ref. [34,35]. 



11 



When a two point function is attached to a vertex, the corresponding time has to be inte- 
grated over. A closed loop corresponds with an integral over spatial momentum J d^p/ (27r)^. 
The tree level contribution to the equal time two point function 

j d^x e-'^-^'itpir, x)^(t, 0)), (60) 

is given by 

Fik,T,T). (61) 

There is no contribution with the two point function because that vanishes for equal 
times. 

It has been shown in Ref. [33], that on large scales, in the classical approximation, we 
can ignore the vertices in Eq. fl59|) with three dashed lines. 



4.2 Non-perturbative tri-spectrum 

In the classical limit where we can ignore the vertex with three dashed lines, the tri-spectrum 
will now simply be given by 

3A 

{^pk,.--i^k,) = + + (62) 

n dr]' 

X J -ij^^F{v,v'',h)F{r],r]';k2)F{ri,T]';k3)G^{T],r]';h) +perm, 

where perm symbolize all permutations of 1, 2, 3, 4. 
Using 

F{x,y) = Q{G-+{x,y)) (63) 
G^{x,y) = 9{xo-yo)2^{G-+{x,y)), (64) 

one can find the leading order contributions to the Green functions in the small momenta 
limit, by expanding the Hankel funtions in Eq. psjl in the small argument limit. In this way 
one obtains the infrared limit of the propagators 

F{k; n, r2) ^ ^{-knYi-kr^Y, (65) 

and 

G^'ik; n, T2) ^ e{n - T2)^ [i-knYi-kr^f'' - {-k^Yi-knY-')] . (66) 

For convenience we have defined 6 = m\^/'iH^. With these expressions, the non-perturbative 
solution for the tri-spetrum in the large wavelength limit becomes 
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where in the last step we took the hmit of small 6. One may worry if we should take 
the limit 77 — > before taking the hmit of small S. However, in the integral above the 
non- vanishing contribution comes from the super-horizon part, which implies that in the 
expression above kr] = —{k/H)exp(—Ht) is exponentially suppressed by the number of e- 
foldings left of inflation after the mode k has left the horizon. For physically interesting 
modes the number of e-foldings left of inflation after horizon crossing is less than about 60. 
Thus, taking immediately 77 — > is not entirely physical, and the limit of small 5 in the 
equation above is a good approximation when 5 ^ 1/60. 

Inserting the value 5 — m^p/Si?^ = •\/3A/127r finally yields 

{A, . . . i>..) « '-^S'^Hk, + . . . + S,)H'Mi^±|^, (68) 

SO effectively the tri-spectrum experience a non-perturbative enhancement by a factor 127r/ a/A 
when A is small. Notice also that the higher-loop corrections to the trispectrum are sup- 
pressed in the IR by the fact that the propagator F{k,Ti,T2) scales like k"^^'^ in the IR. 



5 Non-perturbative corrections to the curvature per- 
turbations 

Now, let us as a final example, discuss an example similar to the one discussed by Weinberg in 
Ref. [12], where the interaction between the curvature perturbation ( and a set of test scalar 
fields (pn was considered. If we assume that the test scalar fields arc only gravitationally 
coupled to the inflaton, the part of the Lagrangian containing the test scalar fields and their 
gravitational interactions is, similarly to in Ref. [12], given by (with SttG = 1) 



n n 

(69) 

where the spatial part of the mertic in this gauge is given by 

g'ij = a^e^^[exp7]jj, 7^ = 0, 9^7^^ = . (70) 

where 7ij(x, t) is a gravitational wave amplitude, and C(X)^) is the curvature perturbation. 
The other components of the metric are given by the lapse and the shift functions, A/", and 
jV*, in the usual ADM approach 

^00 = -hf" + 9io^'M' , c/io = 9i3>^' ■ (71) 

Now we can expand A/" in a sequence in perturbation theory J\f = 1 + 0;'^^) + 0;^^^ -|-. . . . Taking 
the potential of the 0{N) invariant test scalar \^(0) — (A/8A")(0 • 0)^, then from expanding 
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the exponential, we would naively conclude that there should be an interaction term of the 
type 



9 



3^2/ 



(72) 



although one can not trust the slow-roll order of that term. Since the Lagrangian for the 
curvature perturbation should vanish in pure de Sitter limit, the Lagrangian should be 
suppressed by some power of the slow- roll parameter [36,37]. To find the right slow- roll 
order of the interaction term, it is convenient to use instead the variable C„ [36,37], which is 
defined as e^n = where Q is the infiaton fluctuation in the uniform curvature gauge and 
e is the first slow-roll parameter. Writing the Lagrangian in the uniform curvature gauge to 
second order in one finds a slow-roll suppressed contribution of the form0 

1 



32N 



(73) 



From the interaction above, we can calculate the contribution to the two-point function of 
curvature perturbations 



(CnCn) 



(27r 
1 



x^e'XNG>ix,x)G>ix,x). 



Using 



and 



we obtain 



Gl{k;,7^o,v)=^-Uk{r]o)U:{ii) 



iH 



(CnCn) 



kV2k 



[1 + ikifjc 



-ikrj 



(74) 

(75) 
(76) 



MiVoM'iv)] G>{x,x)G>{x,x) 



(2vr) 



12^ 



y + ^Ci(-2A;r^o) + 



G^{x,x)Gl{x,x) 



(77) 



If we insert the asymptotic value G^{x,x) = G^^{x,x) ~ H'^/2tt\^, which we calculated 
in section 4, then we obtain (remember we are taking SttG = 1) 



(CnCn) 



AttNH^ 

6F 
87r2 



20 



3A;3 



eNH'n 



- -Ci(-2A;r/o) + . 
y + ^Ci(-2A;r/o) 



(78) 



■l-The term proportional to a"^^V{(^^ cancels out after applying the background equations. The leading 
surviving contribution is proportional to (a^^'')^F(0), where a^^"* = (Xj'^e^Q^ was given in Ref. [36] (while 
a\ was computed in Ref. [10,38]). 
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where is the tree-level power spectrum of curvature perturbations and the logarithmic 
running is because C„ is not conserved. However, we can evaluate the expression at horizon 
crossing when Ci{—kr)o) ~ 1 and make a straightforward transformation into the conserved 
Because wc are considering initially a two-loop correction, wc see that there is an suppression 
of an additional which typically will make the contribution small. Nevertheless this 
simple analytic solution illustrates the behavior of non-perturbative effects, and, like in the 
previous section, we find a parametric enhancement in the coupling constant A. One may 
then speculate about the non-perturbative behavior of other more realistic models with one- 
loop contributions, which will scale with one power of less and where we know that the 
perturbativc one-loop result can already be significant. However, a detailed study of such 
more complicated cases are beyond the scope of the present paper and is left for future work. 

6 Summary and conclusions 

In this paper we have analyzed the properties of the cosmological perturbations during a de 
Sitter stage of a set of scalar fields subject to an 0(N) symmetry with quartic interaction 
and in the limit of large N. In such a regime the theory can be solved exactly keeping control 
of the IR cumulative effects which appear in perturbation theory and are proportional to the 
total number of e-folds of the infiationary stage. One finds that the correlation functions are 
IR finite. In particular, we have computed the resulting trispectrum, which is the interesting 
quantity signalling the presence of NG's. We have shown to parametrically differ from the 
tree-level value. Its amplitude is suppressed only by the square root of the quartic coupling, 
thus manifesting the presence of non-pertrubative effects. In our study we have limited 
ourselves to the computation of the scalar field correlation functions; we have not addressed 
the issue of relating these correlations to those of the truly gauge-invariant (iso-) curvature 
perturbations. However, this step is automatic through the SN formahsm [39-43] once the 
time-evolution of the universe is known. Finally, we have computed the non-perturbative 
corrections to the two-point correlator of the comoving perturbation coupled to a seto of N 
fields, i.e. the toy model proposed by Weinberg [13] to show the presence of one-loop IR 
divergences. We have shown that the non-perturbative corrections are again parametrically 
different from the one-loop result. Our findings may be considered as a first step towards 
the understanding of the IR cumulative effects which plague at the loop-level cosmological 
scalar perturbations when the theory is self-interacting and a superlarge box is assumed. 
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